Macroscopic Degeneracy and Emergent Frustration in a Honeycomb Lattice Magnet 
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Using a hybrid method based on fermionic diagonalization and classical Monte Carlo, we investi- 
gate the interplay between itinerant and localized spins, with competing double- and super-exchange 
interactions, on a honeycomb lattice. For moderate superexchange, a geometrically frustrated tri- 
angular lattice of hexagons forms spontaneously. For slightly larger superexchange a dimerized 
groundstate is stable that has macroscopic degeneracy. The presence of these states on a non- 
frustrated honeycomb lattice highlights novel phenomena in this itinerant electron system: emergent 
geometrical frustration and degeneracy related to a symmetry intermediate between local an global. 

PACS numbers: 71.10.-w , 75.10.-b , 71.27.-|-a , 71.30.-|-h 



The Kondo lattice model (KLM) is probably the most 
celebrated starting point for the investigation of the in- 
terplay between localized spins and itinerant electrons 
[l| . ft provides the canonical explanation for the Kondo 
effect and for the heavy-fermion behaviour observed in 
many materials '2] . Motivated by the search for topolog- 
ically non-trivial states of matter, several groups have 
recently studied the itinerant KLM on frustrated lat- 
tices, such as the triangular or the pyrochlore one, and 
have shown that due to the strong geometrical frustration 
scalar-chiral types of magnetic ordering emerge 0-11] . 

The physics of the KLM on non- frustrated lattices, 
such as the square and cubic one, has been studied 
extensively. In particular the limit of strong coupling 
and large localized moments, where the KLM goes over 
into the double-exchange (DE) model, is directly rele- 
vant to the colossal magnetoresistance effect in perovskite 
manganites In such cases, the competition be- 

tween DE and antiferromagnetic (AFM) superexchange 
can lead to canted spin states or phase separation Q . Al- 
though the honeycomb lattice is also bi-partite, it has the 
smallest possible coordination number for proper 2D lat- 
tices. That the honeycomb lattice can support physical 
phenomena fundamentally different from square lattices, 
is illustrated by recent Quantum Monte Carlo calcula- 
tions , which identify a novel spin-liquid phase for the 
Hubbard model on the honeycomb lattice, a finding sup- 
ported by analytical studies [lol - [l^ . 

In this Letter, we investigate the consequences of the 
competition between AFM superexchange and ferromag- 
netic (FM) DE on the honeycomb lattice. We find that 
two exotic ground states exist between the trivial, fully 
FM and AFM phases. In the first, nearer to the FM 
state, the spins self-organize into FM hexagons that are 
coupled antiferromagnetically. Since the hexagonal rings 
form a frustrated triangular lattice, their order is remi- 
niscent of the Yafet-Kittel state [13[. The competition 
between isotropic magnetic interactions thus causes geo- 
metric frustration to emerge in a non-frustrated lattice. 



For slightly stronger AFM interactions, wc find the 
exact groundstate to consist of independent FM dimers 
containing one electron each. Apart from the require- 
ment that the alignment of adjacent dimers be AFM, 
they are independent. The groundstate of this TV- 
spin system therefore has a high degeneracy oc . 
While the macroscopic degeneracy oc e"^ in (spin) ice 
is caused by the local symmetry of the frustrated tetra- 



hedra [lj,ll5|: our exponent indicates the presence of 
an 'intermediate' symmetry - a symmetry between local 
and global fid\ . It is remarkable that this highly degener- 
ate groundstate manifold arises as an emergent effect in 
a Hamiltonian that itself does not have such a symmetry. 

In many materials, the essence of the electronic struc- 
ture is captured by interacting spins and electrons on a 
honeycomb lattice. Interactions between impurity mag- 
netic moments on the honeycomb lattice of graphene 
have been intensely studied in a Ruderman-Kittel- 
Kasuya-Yosida (RKKY) framework 17] and using the 
KLM fis'l . Going beyond RKKY is even more important 
in transition metal oxides, e.g., Bi3Mn40i2(N03) [l^ or 
Li2Mn03[23], with Mn ions on a honeycomb lattice. 

The Hamiltonian corresponding to the one-band DE 
model in the presence of competing AFM superexchange 
interactions on a honeycomb lattice is 



iJ - - 2^ t,j 



+ H.c.) + JafY^S, 

{ij) 



S„ (1) 



where V'j and "0^ are the fermionic creation and anni- 
hilation operators, respectively. In accordance with the 
DE scheme these fermions have their spin aligned with 
the on-site spins S^. The on-site core spins are treated 
as classical spins with |Si| — 1 and thus can be speci- 
fied by their polar and azimuthal angles (0^,0^). Both 
sums are over nearest neighbors. Due to the alignment 
of electron spin to the core spins, the hopping ampli- 
tude depends on the direction of the core spins, tij = 
io[cos(6li/2)cos(6'j/2) -Hsin(6',/2)sin(6'j/2)e-'('^'-'^')] 0. 
The strength of the AFM super-exchange is given by Jaf 



2 




b) 



A-sublattice oB-sublattice 



FIG. 1. Schematic view of (a) the honeycomb lattice and (b) 
the brick-wall lattice having the same topology. 



and all energies are in units of the hopping amplitude to ■ 
To guarantee an unbiased search for groundstate can- 
didates, we employ a well-established hybrid method of 
exact diagonalization (ED) for the bilinear fermionic part 
of the Hamiltonian and Monte Carlo (MC) for the clas- 
sical spins [01 . Each MC configuration is defined by a 
given core spin texture and Markov chains are generated 
by diagonalizing the fermionic problem for each configu- 
ration update. We also make use of the travelling cluster 
approximation (TCA), which has proven its validity and 
success in earlier studies on a similar class of models 8] , 
to go to larger lattice sizes. We report here results based 
on calculations on an = 12^ honeycomb lattice, using 
a cluster of size Nc = 6^ . In the MC routine we use ~ 10^ 
steps for equilibration and the same number of steps for 
thermal averaging. We focus on the case of a half-filled 
band, which refers to 1/2 an electron per site, equivalent 
to quarter-filling in the spinful problem. For selected pa- 
rameter values, the MC procedure was further refined by 
an optimization routine that diminishes thermal fluctua- 
tions (21]. 

To identify the magnetically ordered states, we cal- 
culate the spin structure factor S'(q) = J^i j (^i 



where 



) is a thermal average and 
is the position space vector of site i. For a clear un- 
derstanding of the real-space structure of the magnetic 
states it is helpful to look at S{q) on a square geom- 
etry [see Fig. lUb)]. A specific long-range ordering is 
expressed as the point in the Brillouin zone where the 
structure factor shows a peak. To analyze the electronic 
properties we compute the density of states (DOS) as 
D{uj) = {jjJ2k^i^ ~ ^k)) and approximate the delta- 
function by a Lorentzian with broadening 7. 

In the absence of super-exchange interaction (Jaf = 
0), the spins order ferromagnetically, as expected from 
the DE mechanism. The fermionic problem is then equiv- 
alent to non-interacting spinless electrons on a honey- 
comb lattice, giving rise to a dispersion and DOS that 
is well-known from graphene [see Fig. HJa)]. Introduc- 
ing a small Jaf still leads to a FM ground state. At 
Jaf « 0.14, the FM state becomes unstable and gives 
way to a state with S{q) peaked at § (tt, 0) (and the points 
related to it by symmetry), and with the peculiar four- 
peak DOS shown in Fig.lUJb). Real-space snapshots show 



that a superlattice formed of hexagons emerges at low 
temperatures T, as depicted in Fig. |3l[a) . This result was 
corroborated by zero-temperature optimization of the 
spin pattern. Spins within one hexagon are almost FM, 
the allowed energies for electrons moving on a six-site 
ring are — 2tocos0 = —2to and — 2toCOS7r/3 = —to, with 
twice as many states at —to, which gives precisely the 
DOS seen in Fig. [2][b). Coupling between the hexagons 
is AF, but since they occupy a frustrated triangular lat- 
tice, see Fig. EJa), perfect AFM order is not possible. 
The hexagons instead are at an angle of ~ 27r/3, cor- 
responding to the Yafet-Kittel state [13| well known for 
the triangular lattice, leading to the signals at |(7r,0) in 
S{q). Thus a geometrically frustrated triangular lattice 
emerges spontaneously from isotropic, competing inter- 
actions on the non-frustrated honeycomb lattice. 

For 0.18 < Jaf ^ 0.25, we find a state consisting of 
classical dimers. The dimers each consist of two spins 
aligned in parallel, they cover the lattice in such a way 
that the neighboring dimers are anti-parallel with respect 
to each other. In Fig. [S^b) and |3l[c) we show two possi- 
ble dimer configurations. In this spin texture, the elec- 
tron kinetic energy reduces to that of uncoupled two-level 
problems, having only two eigenenergies ztto. The DOS 
is therefore given by £)(aj) = S{(jj — to)/2 + S{uj + to)/'2, in 
excellent agreement with MC calculations [see Fig. [D^c)]. 
The dimer state can be understood as a trade-off between 
the FM ordering and the AFM ordering: the electrons are 
allowed to populate all the —to levels (which is more fa- 
vorable compared to AFM) and the spins are anti-parallel 
with respect to two of their nearest neighbors (which is 
more favorable compared to FM). 

Interestingly, the dimer ground state of this quantum 
system has a macroscopic degeneracy, i.e., there is a 
macroscopically large number of ways to cover the lat- 
tice by dimers such that the neighboring dimers are anti- 
parallel. One way to see the degeneracy is to start cov- 
ering lattice rows in Fig. [ifb) by dimers. It is easy to 
see that having fixed the dimer pattern in the 1st row, 
there are two independent ways of covering each subse- 
quent row, giving 2^^^ states for a A^-site lattice. The 
fact that there is thus no long-range order along the y 
direction of the brick- wall is reflected in S'(q), which be- 
comes flnite along lines in momentum space, as in com- 
pass models 24|427l|. In the 2D compass model, differ- 
ent degenerate configurations can be reached by flipping 
a row of spins. The corresponding \/N operators com- 
mute with the compass Hamiltonian and thus define an 
intermediate symmetry, i.e., between a local, gauge-like 
(oc e^) symmetry and global one (independent of A^) [l6| . 
The magnetic order parameter that obeys the interme- 
diate symmetry is consequently of nematic type. In the 
dimer state, the minimal symmetry operations involves 
translation of all spins in two adjacent zig-zag rows by one 
lattice spacing, cry 1— (Tij+i [cTy is the spin at site (i, j)]. 
An example for two dimer configurations connected by 
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FIG. 2. (Color online) (a)-(d) DOS at low, intermediate 
and high temperatures for different values of Jaf (7 ~ 0.04). 
In (a) the T = 0.001 curve shows the DOS of free fermions 
on a honeycomb lattice in the thermodynamic limit. In (d) 
the T = 0.001 curve represents a gapped insulating phase, 
the seemingly finite DOS at Ep being a broadening effect. 
Except for the FM phase, all ground states are gapped. 



such an operation is given in Figs. [Sfb) and[31Jc), where 
the second and third rows were shifted. However, this op- 
erator does not commute with the Hamiltonian Eq. ([l}, 
and the intermediate symmetry is thus rather a prop- 
erty that emerges in the system's ground state, similar 
to the case of striped phases at fractional filling in the 
regime of narrow bandwidth and small Jahn- Teller cou- 
pling in a model used for manganites [28| . This inter- 
mediate macroscopic degeneracy should lead to a large 
specific heat at low temperature. 

For strong super-exchange coupling, there is a continu- 
ous way in which the dimer state can approach the AFM 
ordered state, captured by a canting angle 9 [see Fig. 
[Sjd)], which is the angle between the two spins forming a 
dimer in the pure dimer phase. The spins remain antipar- 
allel to those of the neighboring canted dimers. In this 
way, the two-level dimer systems remain uncoupled. The 
hopping amplitude between the two spins in the dimer is 
renormalized by the DE mechanism to tg cos{9/2). The 
DOS for a canted dimer state is consequently given by 
D{uj) = S{uj-tocos{e/2))/2 + 5{Lu + tocos{9/2))/2, as can 
indeed be observed in the DOS for Jaf = 0.50 shown 
in Fig. md). The canted-dimer groundstate has a gap 
oc cos(6'/2) at the chemical potential, which shrinks as 9 
approaches tt for Jaf —^00. At finite T, the two peaks 
widen and merge due to thermal spin fluctuations, lead- 
ing to a metal with reduced band width, see Fig. EJ^d). 
This canted state retains the macroscopic degeneracy in- 
herent to the AFM dimer state discussed above - also 




FIG. 3. (Color online) (a) Snapshot from MC simulations 
supplemented by optimization routines showing an emergent 
triangular lattice (black circles) formed by FM hexagons at 
Jaf = 0.14. Spins within each hexagon are almost FM, a 
small canting angle between groups of three is illustrated by 
shading. The colored spins illustrate the 27r/3-angle order of 
the Yafet-Kittel state. Schematic view of, (b)-(c) two dimer 
states related by a translational symmetry (see text), and (d) 
a canted dimer state. 



this ordering is therefore of nematic type. 

Our results are in good agreement with elementary 
energy considerations. The energy per site varies as 
3 Jaf/ 2 and —Jaf/ 2 for the FM and the dimer states, 
respectively. This would imply a phase transition at 
Jaf ~ 0.15, the FM state is indeed stable for Jaf ^ 0.14 
and the dimers for Jaf ^ 0.18. In between, the emer- 
gent Yafet-Kittel state, with a more complex energy de- 
pendence, is favorable, see Fig. lUJb). The energy per 
site for the canted dimer state is — (2Jaf — cos{9)Jaf + 
to cos{9/2))/2. By differentiating with respect to the 
canting angle 9, one easily obtains that canting becomes 
favorable for Jaf > 0.25 and that the optimal energy is 
then given by — 3Jaf/2 — ip/ (16 Jaf)- This is reflected in 
the behavior of the ordering temperature for the dimer 
state, which starts decreasing at Jaf = 0.25 [see Fig. 
Ha)]. 

The results are summarized in Fig. 2) In the finite-T 
phase diagram Fig. SJa), phase boundaries for the FM 
and quasi- AFM regions are obtained by determining the 
inflection point in the {M){T) and {M){T) (M denotes 
staggered magnetization) curves. The onset of dimer and 
other phases is determined by tracking the temperature 
dependence of the spin structure factor and the charac- 
teristic features in the DOS. Figure Hl^b) compares the 
ground-state energies of the various phases and perfectly 
agrees with the unbiased numerical data, indicating that 
we have identified the ground states correctly. 

In a full quantum treatment of the spin system ad- 
ditional quantum fiuctuations can affect the stability of 
these ordered phases. Here one anticipates the FM build- 
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mediately raising the question which other models share 
such features and how further residual interactions might 
affect the degeneracy, but pertains in particular to honey- 
comb manganese oxides, which form a promising class of 
materials to realize these novel types of highly frustrated 
states harboring macroscopic degeneracies. 

This research was supported by the Interphase Pro- 
gram of the Dutch Science Foundation NWO/FOM and 
by the Emmy-Noether program of the DFG. 



FIG. 4. (Color online) (a) T-Jaf phase diagram at half-filling 
obtained by MC for a 12 x 12 lattice, (b) The energy of the 
various states: "aim. FM" refers to a spiral with the longest 
wavelength supported by the lattice, converging to FM in the 
thermodynamic limit. Similar finite-size effects are reported 
in doped ID and 2D lattices "Hex" denotes the emergent 
Yafet-Kittel order between hexagons depicted in Fig. EJa), 
the energy was optimized with respect to the canting angle 
within the hexagons. "Dimers" and "C. Dim" are the highly 
degenerate FM and canted dimer states, and "AFM" denotes 
perfectly AFM order. The black crosses are energies obtained 
by unbiased MC and a subsequent energy optimization. 



ing blocks (hexagons and dimers) to be robust as they are 
stabilized by a substantial DE energy, and the FM state 
remains an eigenstate of the hexagon (dimer) for quan- 
tum spins. The Yafet-Kittel ordering between the large 
total spins of the hexagons is expected to be more classi- 
cal, and thus robust, than for S — 1/2, where it is found 
for r — ^ [2^. If one can describe the magnetism of a 
dimer state by an effective NN AFM Heisenberg model, 
then this model remains the same if one performs the 
operation illustrated in Fig. [3l The emergent symme- 
try would thus commute with the effective low-energy 
Hamiltonian so that the corresponding degeneracies are 
preserved. 

We conclude that the isotropic double-exchange model 
with competing super-exchange interactions on the non- 
frustrated honeycomb lattice has an unexpectedly rich 
phase diagram with exotic magnetic phases. In one of 
these, FM rings become the essential building blocks, 
which form a frustrated triangular lattice and are an- 
tiferromagnetically coupled. The stabilization of such 
frustrated spin states on a bipartite honeycomb lattice, 
without explicit frustration, is so far unique and an exam- 
ple of geometrical frustration emerging from competing 
interactions. Another novel phase consists of FM dimers 
ordered antiferromagnetically and has a degeneracy. 
This is reminiscent of compass models, but in the present 
case the corresponding symmetry is not a property of the 
Hamiltonian given a priori, but rather a property that 
emerges in the systems ground state 2J-l28j . These phe- 
nomena are not only relevant in a theoretical context, im- 
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